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Abstract 

In this paper we study the chiral symmetry breaking in the hard wall AdS/QCD model. We solve 
OO , the equations of motion up to the second order at large momentum and obtain the first few terms in 

■ the expansion of the "left-right" correlator, which is the measure of symmetry breaking. We study the 

' dependence on AdS radius to get the result as the series in t'Hooft constant. 

H 

i3 ' 1 Introduction 

In the last few years a great attention was paid to the so-called phenomenological AdS/QCD theories. The 
essence of these models is to use the AdS/CFT correspondence [IJ to describe QCD in large Nc limit via 
its 5-dimensional dual theory. The exact structure of this 5D theory, describing all specific features of QCD 
is not clear, but some simple models have been proposed [2] [3] [1] [S] [S] which already give promising 
results. 

In this paper we study the simplest of these settings, the so-called hard wall AdS/QCD model (see 
' for example [2] [1], first proposed in [8J) . Our goal is to find solutions to the equations of motion in 5D 
(~| I theory perturbatively in inverse powers of momenta and analyze two-point correlators at large and their 
dependence on AdS curvature radius. Classical solutions play an important role in AdS/QCD, because, 
I according to AdS/CFT correspondence, correlation functions in QCD can be evaluated as variations of 5D 
^ ' effective action on the classical trajectories with respect to the boundary values of the corresponding fields. 

vn . 

Z^^[Jl{x), J2{x), ...] = exp {Slff[Vi{z,x) classic, ^2(2, x)classic, •■•]) \Vi(Q,x) = .h(x) 

CN ■ 

\l \ This rule allows us to compute 2-point functions for vector, axial-vector and pseudoscalar currents in 

QCD, fix the free parameters of the model and also study the so-called "left-right" correlator {{LB) = 
\ {^y) ~ (^^)) up to the term, quadratic in condensate. The dependence on QCD coupling constant is 

reconstructed in the result, via the ADS/CFT correspondence of Yang-Mills coupling constant and AdS 
curvature radius (see for review |16j). 



Oh: 



00 
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We obtain the expression for the "left-right" correlator at the strong coupling regime. The first term, 
contrary to the sum rules result |17| , which is linear in A', is proportional to A'". This is not strange, because 
calculations via AdS/QCD imply the large coupling constant, while sum rules work well at weak coupling. 

This paper is organized as follows. In Section 2 we give a review of AdS/QCD hard wall model and 
derive action for axial, vector and pseudoscalar fields and corresponding equations of motion. In Section 3 
we solve these equations in the limit of large momenta by means of the Green function, derived in Appendix. 
Section 4 is dedicated to the evaluation of the correlation functions via the AdS/CFT recipe and matching 
of parameters in result with QCD ones. Conclusion is given in the Section 5. 
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2 The description of the model 



Here we consider the simplest holographic model of low energy QCD, proposed in [2] , [4] , [6] (see also [131 
UM [T5]). the so-called "Hard wall AdS/QCD model". In what follows we will work with conventions and 
notations, used in [5]. 

In the AdS/CFT prescription, the fields in 5-dimensional space are dual to operators in 4D, and the global 
flavor symmetry of the 4D field theory corresponds to the gauge symmetry in its 5D dual. So we will study 
4D QCD with SU{2)lX SU{2)r global symmetry via the gauge theory in AdS with SU{2)l x SU{2)r gauge 
group. In this model only the fields, dual to QCD operators with the lowest dimensions, are considered. 



2.1 The metric and the fields 

We have the SU{2)l x SU{2)ji gauge field theory in AdS^ space with the metric: 

7->2 

ds^ = ^{-dz^ + dx^dx,,) (1) 

where R is the AdS curvature radius, cut at z coordinate: < z < Zm 

Later, we will denote 5-dimensional indices with capital Latin letters {A,B.. G {0,1,2,3,2;}), and 4D 
indices with Greek letters (/i, i^, .. S {0, 1, 2, 3}). We will use the metric tensors, defined as follows: 



9AB 



diag ( - 



p2 n2 p2 7j2 td2 
it it it It It 

\ z^ ' z^ ^ ' ' 
VAB = diag{-l, 1, -1, -1, -1) 



g = det{g_ 



AB) = 



The theory includes left- and right-handed gauge vector fields of SUl{2) x SUr{2) {Al and Ar, 
respectively) and bifundamental scalar Xap. According to AdS/CFT 5D fields correspond to operators 
in QCD: 

with the boundary conditions imposed at z — Zm'- 

a.y(z™) = o ; a,A(z„) = 0. 



The action is: 



where 



d'x^Tr !^A'i\DX\' + ^\X\') - ^^{FI + F'^) \ (2) 
DbX = ObX - iAlbX + iXArb 



Fbd = Ob Ad - do As - i[Ab,Ad], 
and we introduce the normalization constant A of field X. 

From lagrangian for X in absence of A and V we get the equation of motion: 

^3x = ^d,d^x-d,^d,x 

z° z'^ z'^ 
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which has the sohition: 

Xo(z) = ^Afz + 

According to AdS/CFT [IB], we argue, that M corresponds to quark mass matrix, i.e. the source of operator 
q'^q^ and E to condensates, i.e. vacuum expectation value of q^q^. We can make M to be quark masses 
exactly by appropriate definition of normalization A, but the relation between E and condensates is to be 
ascertained. In further discussion we set E = crl, M — ml, assuming the equality of quark masses. 

Xo{z) ^ ^v(z)l, v{z) = mz + az^ (3) 

2.2 The action for vector and axial vector fields 

It is convenient to rewrite the action in terms of vector and axial vector fields: 

V^{Al+ Ai^)/2 
A^{Al-Ar)/2. 

To get the effective action for small fluctuations of fields, we will derive it up to quadratic order. 

Fl = F,,j,sFt^ = F,.,,sF,M^V^'')i^Vn 

Fl.abFl.cd = {OaVb + OaAb - ObVa ~ ObAa) ■ {dcVo + dcAo ~ doVc ~ doAc) - 

= Fv.abFv,cd + Fa^abFa.cd + Fv^abFa,cd + Fa^abFv.cd 
Fr.abFr^cd — Fv.abFv^cd + Fa^abFa.cd — Fv^abFa,cd — Fa^abFv.cd 
As can be seen, the cross-terms cancel, so the gauge part of ^ takes form: 

Saa,vv = I d^x^ ^-^^ i^v'^'^nFl^BFlcD + F^,ABF^^cD)Tr{t^t''} = 

= I d'x A j + F^F^) (4) 

where we introduced the notation: •q^^rj^^ F'^ abFa cd — FaFa- 

Note, that in F\ and Fy we have the same terms with odd number of A, arising with opposite signs. 
So they will cancel each other and there will be no terms in result, containing odd number of axial fields A 
(for example AVV). 

Now we consider the part of the action ([2]), describing the interaction of X with the gauge fields and 
decompose X as follows: 

^ X^^exp{i2Tr''{eyl^), X^'' = (5""^ ■\{mz + crz^) = (J"'' • ^v(z). 

Expanding X to the first order in tt = 7r"<°: X = Aro(l + 2i7r) we get 

DaX = Xo{2idATT - iLa + 2Latt + iRa - 2ttRa) 
DaXDbX = v{zf{dTT - A)A{dTT ~ A)b 



therefore 



S^r^t^ J d''x^Tr{A^\DX\^} (5) 
= J d'x^i^T^^^)AMz)'{dn'' - A'')A{dn'' - A'')BTr{et'} := 

d'x ^-^{dn-Af. (6) 
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2.3 Equations of motion 

We have obtained the action in terms of V, A and tt. Now, let us derive the equations of motion, imposing 
the gauge 

A, = V, = (7) 
d^V^ = 
Af, = A^^, + df,(j) df,Aj_f, = 0. 

For Va we get (g]): 

z z 

z z 

and after the Fourier transform V{q, z) — J d'^x e''^y(x, z) it takes form 

d^d,V^-iq,z) + ^V^^{q,z)^0. (8) 
The action of axial sector in terms of A±^, A^, <j) and vr is (H]),®: 

^/^BM^ [-(S.TT - A.r + id,- - d,<t>f + {A^,r] (9) 
and gives the following equations of motion (the gauge = is used): 



= for A^^, (10) 



± 



z z z-^ 

Qa.r) + ^'g^A^^^ (tt" - 0°) ^ fora^<^, (11) 

-q^d,cj,-+^M^d,T:- = Q forA,, (12) 

2 2 

a.^9.^ + g2^(^-0) -0 forvr. (13) 
One can see, that pip follows from (fT^ and (fT5|) . so we can use the last two to find solutions for and tt. 

3 Solution to the equations of motion 

According to AdS/CFT correspondence, the generating functional of correlation functions in 4D conformal 
theory equals to the effective action of its dual 5D theory computed on the classical trajectories [T] [16j. 
Correlation functions can be obtained by the variation of the 5D action with respect to the boundary values 
of fields. For example, vector current correlator is: 

S S 

{Jv{qi)Jv{q2)) = . . . 5'(Vb)|vo^o- 

oVo[qi) oVo(q2) 

It is convenient to introduce functions v{q, z),a{q, z) and (j){q,z), such that: 

VA{q, z) = VoA{q)v{q, z) ; w(q, z%=^ = 1 

Aa{sI,z) = AaA{.q)a{.q,z) ; a{q, z%_=^ = \ 

0(9,z)U=. -(/<o(g) (14) 
and use the euclidian momentum = —q^. 
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3.1 Solution for V 

The equation for the function v ([5|) P^ is: 



Oz-UzV V — 



z z 

which can be reduced to: 



«" + — +■« -0^-^ -0, u=-. 
z \ J z 

This is the modified Bessel equation with A = 1. Its sohition is: 

u{z) = const{AIi{Qz) + BKi{Qz)) 

where /i and Ki - are modified Bessel functions of the first and second kind respectively. Constants in this 
solution can be fixed, using the boundary conditions on v{z) at z = e (1141) and z — z„i ([7]), 

A = Ko{QZyn) B = -Io{qZm) 

const = ^ 
B 

and we obtain the solution for v. 

v{Q,z) = - wi . Qz[Ko{Qz^)h{Qz) - Io{Qz^)Ki{Qz)] = ^x{Ah{x) + BKi{x)). (15) 

3.2 Solution for A 

The equation for A±_ (fTU]) differs from the equation for ([5]) by the interaction term. The solution can't 
be found analytically, but we can study its asymptotic at large euclidian momentum Qz^ ^ 1 in the chiral 
limit w? 0. We will also keep the ma term, to study the condensate dependence of the solution. 

Oz -Oza a ^ a 

\z J z z-^ 

--dza + dla - Q^a ^ R^glK^a^z'^a + R^glK^{2mcr)z^a 
Changing the variable Qz ^ a;, we get at large Q the inhomogeneous Bessel equation (A = 1) for function 



Q ^ (J Q i\ ('2.(7 Tfi\ Ft 

xa{x) (for convenience we denote small parameters A = qs — - and /i = — — q4 ' — ): 



dxa + d'^a — a = Xx'^a + j^ix^a (16) 

It can be solved by means of the Green function, derived in [Tj, which for the euclidian momentum is 
presented in Appendix. 

We will compute first and second order corrections due to A-term, because they don't vanish in the chiral 
limit and derive ^-term, related to quark mass. 

3.2.1 First order correction 

The first order correction to the solution to the homogeneous equation can be computed, using the Green 
function ([33]) . as the integral (xq — Qe, Xm = Qzm)' 



(1) = fdx'Xx'^a^'Hx')^^^^ 
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where a^^^^ - is the solution to the homogeneous equation (fT5|) . found in the previous Section. 
Inserting the Green function (l33|) . we have: 



= -^^'f^+^^^(-)^ £ dx'Xx-aio^ix') . x'iChi^') + DK,ix')] + 
^xjCh^+DK^ £"dx'Ax'3a(°)(x')-x'[M(x')+i?i^i (.')]. (17) 

We are mainly interested in the behavior of a{z) at the vicinity of the boundary z ^ (consequently x — > 0), 
hence the major contribution is due to the second integral: 

^ x[Ch{x)+DK,{x)] r- ^^,^^^,3^(0)^^,) . 5„(0)(^,). 
^ ' AD - BC Jq ^ 

In the limit xq = qe and a:; — > we get: 

2 Xq Xq 

1 . . 



therefore 



(18) 



In the limit '"^ (x) (fTC]) takes form: 

VV 2a;,„ V 27ra;„^ / 

and the integral equals to: 

00 00 

J dx'x'3[a(0)(a;')]' = J dx' x'^K^ix')]^ = | 

Consequently we have the first order correction to the solution (near the boundary) : 

a^'HQ,^)--liQ^r'-^^^- (19) 

3.2.2 Second order correction 

To compute the second order correction we have to evaluate integral: 

a(2) = 7dx'Xx'\<^'\x')^^^. 
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Inserting the Green function (l33|) . we get: 



x[Ah{x) + BKi{x)] 
AD - BC 



dx'\x'^a^^^x'[CIi(x') + DKi{x')]- 



x[Ch{x)^DK,{x)] I ^,,y^^a^ii)^,^Aj^^^>~^^BK,{x')]. 
AD — BC I 



Similarly to the first correction calculation, we are mainly interested in the behavior of a{z) near the boundary 
z — > (and a; — s- 0), so the major input is due to the second integral. We see, that to compute this integral 
we have to know the form of the first correction at any x (not only near the boundary, where we've found 
it ([TO]) ). Therefore we can try to find the form of the function a'^^^(x) (and its order on Q and Zm) near the 
boundary and evaluate the coefficient numerically. 

In the limits Xm ~> oo and xq — + 0, the coefficients A,B,C,D behave as: 



so a(^' (dH) takes form (fTH]) : 



^1 
B\ 



Ko{Xjn)\x^->c 
-Io{Xm)\xm^ 



2r 







V^nXr. 



CUo^o = -^i(a;o)Uo-^o = — 

Xq 



a'-^\x) ^ xKiix) J dx'\x'^Ki{x')h{x') + xIi{x) / dx' \x'^ Ki{x' ) Ki{x' ) 





and we can rewrite a'^'^\x) as: 



where a is the numerical coefficient: 



a^'^{x) 



(20) 



a = j dx'x"^Ki{x') 





Ki{x') j dx"x"''Ki{x")h{x") + h{x') / dx"x"''Ki{x")Ki{x") 



i 124. 



3.2.3 Quark mass correction 

The /x-term correction can be easily computed similarly the subsection 3.2.1 via the integral: 



(1) /" J ' '2 (0)/ /^ C!{x, x ) 
\J = dx fix '(a; ) 
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and after inserting the Green function and taking the appropriate limit: 

oc 



dx'x'[a^°Hx')]^ = / dx'x'^[Ki{x')]^ = I 



The integral equals to: 



and we get: 

^^M--4 '''''^'>"' . (21) 

Finally, from p5ll9l20l2ip we obtain the solution to the equation of motion for a(z) up to the second 
order in A at the vicinity of the boundary (z 0) at large Q^: 

«(^) = wi Oz[KoiQzra)hiQz) - UQz,^)Ki{Qz)]- 
Io[Qz„i) 

3.3 Solution for and n 

Now we will compute to the leading order solutions for the pseudoscalar fields in the theory. These are 
solutions to the equations (fT^ . p^ with fixed boundary value of ^ at z = e JH]). Differentiating and 
substituting dzcj) from IT^ we get the equation: 



dl^d^TT - [dz^] -ird^^d^n - Q'^—d^n ^ — g^-^ ^ q. 

\ z'^ J z-^ z-^ z° 

We need to solve it near the boundary, so we can substitute for v{z) its asymptotic v{z) = mz\z^a Q and 
after change x = Qz it takes form: 

d^-d^TT + -d^-d^TT 9a; TT d^TT = 0. 

X XX X Q'^ X 

In the large limit we can neglect the last term and obtain the modified Bessel equation with A = for 
the function ^dxir. Hence the solution for 7r(z) is: 

^(z) = A'Qzh{Qz) + B'QzKi{Qz). 

Then, using (fT^ we immediately obtain the solution for (f>: 

The boundary condition on at z = e (fT4| fixes the constant B': 

r\2 



(z) = --^^^ Qz[j^h{Qz) + B'K,{Qz)]. 



^' " 2P2A2 2 '^0(g) 



and we finally get solutions at z ^ 0: 

<^(z)|2=e = 0o(g) 

7r(z)|3=, = -~0o(g) 2P2A2 2 



5.<^(z) 



z 

a^7r(z) 



^o(g) .2iL. ?^»(Q^^^)- (23) 



55 
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4 Computation of correlators 



It was already mentioned in the previous Section, that two-point correlation function can be obtained by 
the variation of the effective action with respect to the boundary values of fields. 

4.1 {Jv, Jv) correlator 

For the beginning let us evaluate the variation of action on the classical solution: 



SS = I d^x^SV + t^SOmV = / d^xdz 
oV oomV 



5L 6L 



5V- 



6d,V 

Hence on the classical solution the variation of action reduces to the boundary term and we get for V^: 



5Sv 



d^x 



R 



9i 



5zV^((72, z) 



Taking twice the variation with respect to the boundary value Jq) (fT4|) we obtain the current correlator: 



v{q, z 



dyv{q,z) 



^"''(<ZM'Z'v-g^9M.)ny(<z^) (24) 



and for nv/(g ) we have, using the solution for v (fT5 



dzv{Q,z) 



Qz{KoiQzrn)IiiQz) - k{Qz^)Ki{Qz, 



Z \lo{QZm) 

— g2 (k^{Qz^) - Io{Qzrn)[ln{Qe/2) + j 



lniQ\^) 



Q^oo 



R 



dzv{Q, z) 



R 



(25) 



It was shown in pi, 1] that this result can be compared with the correlator of vector currents in QCD to 
obtain the value of 35. In QCD we have [T7] 



nv(Q2 



"247r2 



hence 55 is set to 



51^ 12^ 
R N, 



(26) 



4.2 {Ja, J a) correlator 

Now we shall evaluate the correlator of the axial currents, associated with A± . It can be evaluated analytically 
only in the limit of large , where we have computed the solution to the equation of motion. 

Using the same procedure as for the vector current we obtain the variation of action on the classical 
trajectories as the boundary term: 



R 



] I d xd qid 

9i 



5Aji_f_,{qi,z 



dzA±^{q2, z) 
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It has the same form as for the vector current, so we have 



Substituting the solution (j22|) . we get 



5.a(o)(Q,z) 



-1 



loiQZm) 



dza{Q,z) 



Q'iKoiQz^) ~ k{Qzra)[ln{Qe/2) + 7] 



Q 



5.(a(i)+a(2)) 



-A - 124A2 + -fi 
5 3 



and finally 



nA(Q') 



JR_ 



2 

2 E>2 A 2 



where A — ^'^^rv, " and /i = ^'''^ ^i^""^ are small parameters in the limit — > 00. 



(27) 



4.3 ( J^, Jtt) correlator 

In this Subsection we will compute the correlator of pseudoscalar currents Jjr = 9759. To obtain this 
correlator, we must find out, what field in our theory is dual to J^r- We have pseudoscalar field (/> = dfj_A^, 
which is connected with axial vector field A^t, dual to the current J a = <il5l^j.<l- Now recall, that: 

If we rewrite this expression in terms of our fields, we find after Fourier transform: 

Q^(j)^2raqJ^ (28) 

and we see, that in our theory the pseudoscalar current is dual to the field 3;^- Hence, to evaluate its 

correlator, we shall vary the action with respect to ^^''^ 

The variation of action of pseudoscalar fields © on classical solutions is: 



R 



Sn—dzTT 



d^xd'^qid^q2e"''''e"''' 



9l 



5(j){qi,z) 



dz(t>{q2,z) 



A^R^m^ 



6n{qi,z) 



dzTT{q2,z) 



hence, the pseudoscalar correlator is given by the second variation with respect to 



{J^{q),JTv{q)) 



4to2 / RQ4 



RQ^ 

2gl glK^A^m^ 



In the chiral limit (m = 0) first term vanishes and we obtain the result: 



R 1 

(J,(<z),J,(g)) -2-2-5^( 



55 9P 



which can be compared to QCD [TZ] : 



{J7r{q),J7r{q)) 



N. 



QCD 



167r" 



This comparison and equation ([26|l allow us to fix A: 

2 _ 8 1 _2Nc 1 



A' 



3 <?|i?2 97r2 R3 



(29) 
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4.4 Evaluation of condensate 



It remained us to fix the relation between a and quark condensate. In QCD one can evaluate condensate as 
the variation of vacuum energy with respect to quark mass. In dual theory we shall variate the action on 
the classical solution. 

5£qcd _ 5 S{Xo) 



5m„ 



5m 



m=0 



The variation of the action ^ on the solution Q is 



6S 



X 



7?3 



f 

/ (f^a; — 5-A^(rn + 3CTZ^)z(5m 
J z-^ 



therefore quark condensate is related to a as: 



{qq) = SR'A'a = 



(30) 



4.5 "Left-Right" correlator 

Now, from (|25p and ([27ll . we can evaluate the quantity Hlb ~ — Hy and after substituting fixed values 
of 35 (P^ . and A (^5]) we get in the hmit of large Q^: 



'l27r2 



64 cr2 64 0-4 

h 124 

15Q6 ^ 9 Q12 



24 aniq 



(31) 



We see, that in this result there is no factors of R, which is by the AdS/CFT [IS] [J correspondence related 
to t'Hooft coupling A' 

i?4 



= asN, = X' = N,gl 



(32) 



47ra'2 ' " 

that means, that in our result all terms are of order 0(A"') in the limit of large A', and we can read out the 
asymptotic forms of coefficients in 

n.« = /(A')^ + .g(A')^ + p(A')^. 



At A' ^ oo our calculation predicts: 



while at weak coupling regime [17j : 



/(A') ^ g{\') ^ p{\') ^ A" 



p(A') ~ A'" /(A') - -4™, ^ A'. 

As for the dependence on iVc, one can see, that provided at large Nc limit {qq) ^ Nc, a ((50)) is of order 
0(iV°), so each term in the result scales as Nc as expected for current correlator. 



5 Conclusion 

From this work we see, that even the simplest AdS/QCD model allows us to compute the current correlators in 
QCD in the subleading approximation. The "left-right" correlator is proportional to (qq) and ruq parameters, 
which violate the chiral symmetry spontaneously and explicitly respectively. In this work we studied the 
dependance of quantities computed on Nc and AdS radius, connected with A'. We reconstruct the power 
expansion of correlators at large and found, that to the leading order all terms are of order 0(A'°). 

The result of this paper is an argument in favor of the validity of the hard wall AdS/QCD model 
in description of chiral symmetry breaking in QCD. It provides additional motivation to search for more 
complicated and suitable models in this branch, moving towards the goal of finding the way to complete 
description of QCD via its holographic dual. 
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A Appendix: Computation of the Green function 

Here we will present the derivation of the Green function for modified Bessel equation ()16p appropriate for 
our problem, which was computed in [71 lllj for Minkovsky momentum. 
Green function is defined as: 

-l + dl~ -dA G{x, x') = xS{x - x'). 



X 

The boundary conditions for G should be chosen as follows ([7]): 

G{xo,x') = dx'G{x,Xm) = 0. 

We will find solutions for x < x' and x > x' and then connect them at x — x' . 
For X > x' let us define variable u = max{x, x'} and at x ^ x' 

'l + dl^^du^ G{u)^0. 

The solution with given boundary conditions is psp : 

G{u) = const ■ u{AIi{u) + BKi(u)) 

A = KQ{xm) ; B = -Ia{x„r) 
For a; < a;' we set v = min{x, x'} and find constants using the boundary condition at w = ccq — > 0: 

G{v) = const ■ v{Ch{v) + DKi{v)) 
b.c: G{xij) = const ■ xo{CIi{xii) + DKi{xa)) ~ 
C = ifi(a;o) ; D = -h{xo) 
Now we can write the Green function in the whole interval in the form: 

G{u,v) = const ■ uv ■ [Ah{u) + B Ki{u)][C Ii{v) + DKi{v)] 

Integrating the equation for G{x, x') by x at the vicinity of x' and assuming the continuity of Green 
function, we get the condition for gluing two solutions. Wherefrom we can define const: 



const ■ x x 



(9,jG'(u, ?;)|u=a;';u=a;' ~ dvG(u^v)\u=x' ■,v=x' — x' 

{AIo{x') + BK^{x')){Ch{x') + DKi{x'))- 

~{Ah{x') + BKi{x')){Ch{:x') + DKa{x')) 

Consider this equation at = — > 0: 

const ■ xq[{AD - BC)Io{xo)Ki{xa) + {BC - AD)Ka{xo)Ii{xQ)\ = 1 
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const ■ xo 



[AD - BC)— + [BC - AD) 



-ln\ — 



2 



= 1 



const 



1 



AD-BC 



Finally we obtain the Green function: 



G{x,x') = 



XX 



AD-BC 



[Ah{u) + BKi{u)][Ch{v) + DKiiv)] 



(33) 



A = Kaixrn); B = -Iq{x„,); C ^ Ki{xo); D = -h{xn) 
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